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INTRODUCTION 


<a 


Recent  trends  in  projectile  design  have  led  to  shapes  with 
greater  length  and  more  slender  ogives.  Unexpected  flight  stability 
problems  have  been  encountered  due  to  decreased  aerodynamic  stability 
of  these  new  shapes.  Clearly,  conventional  aerodynamic  predictive 
capabilities  were  not  adequate.  In  an  effort  to  avoid  these  problems 
in  the  future,  the  Ballistic  Research  Laboratory  has  been  developing 
advanced  numerical  computational  techniques  for  computing  projectile 
aerodynamic  characteristics  to  improve  shell  design  technology.^  — 


Substantial  progress  has  been  made  in  the  past  10  years  in 
the  development  of  aerodynamic  computational  techniques  and  in  the 
availability  of  high  speed  digital  computers.  This  progress  has  made 
it  possible  to  begin  to  use  advanced  finite-difference  computational 
techniques  to  perform  parametric  aerodynamic  studies  for  evaluation  of 
proposed  design  concepts. 


The  use  of  advanced  numerical  computational  techniques  for  a 
parametric  study  is  difficult  to  justify  to  compute  only  static 
aerodynamic  parameters  since  cheaper,  less  complex  techniques  such  as 
Ref.  (1),  (2)  and  (3)  are  available.  However,  if  dynamic  derivatives 
such  as  Magnus  and  pitch  damping  are  considered  important  and  if  viscous 
drag  is  of  interest,  then  the  advanced  computation  techniques  are 
justified  and,  in  fact,  must  be  used.  This  paper  reports  the  initial 
results  of  an  ongoing  research  effort  at  BRL  to  form  an  advanced 
aerodynamic  computation  capability  that  will  provide  the  shell  designer 
with  a  complete  package  of  static  and  dynamic  aerodynamic  coefficients 
for  use  in  design  studies. 
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II.  COMPUTATIONAL  TECHNIQUES 
Scope  of  Effort 

Three  dimensional  finite-difference  flow  field  computational 
techniques  for  inviscid  and  turbulent  viscous  flow  have  been  applied  to 
generate  a  comprehensive  set  of  aerodynamic  coefficients  for 
cone-cylinder  (CC),  tangent-ogive-cylinder  (TOC),  and  secant-ogive- 
cylinder  (SOC)  body  configurations.  The  model  geometries  considered  in 
this  study  are  shown  in  Figure  1.  Body  lengths  up  to  seven  calibers 
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Figure  1.  Model  Geometries 


and  ogive  lengths  of  two,  three, 
and  four  calibers  have  been 
considered.  The  aerodynamic 
coefficients  computed  are 
pitching  moment,  normal  force, 
center  of  pressure,  Magnus  moment, 
Magnus  force,  Magnus  center  of 
pressure,  form  and  viscous  drag, 
roll  damping  and  pitch  damping. 

The  sign  convention  for  the 
pitch  plane  and  Magnus  forces  is 
shown  in  Figure  2.  All 
aerodynamic  coefficients  are 
computed  in  a  conceptually  exact 
manner.  The  only  empirical  input 
is  that  required  for  the  modeling 
of  turbulent  eddy  viscosity. 

The  computations  have 
been  carried  out  for  a  Mach  number 
range  of  1.75  <  M  <  5.  These 
computations  were  all  performed 
for  an  angle  of  attack  of  1°,  a 
nondimens ional  spin  rate  (PD/V) 
of  0.19,  and  for  sea  level 
atmospheric  free-stream  conditions. 
Specific  comparisons  to  wind 
tunnel  data  were  made  for  the 
tunnel  operating  conditions. 


Figure  2.  Magnus  and  Normal  Forces  on  Spinning 
Projectile 
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Coupled  Inviscid-Viscous  Computations 

The  sequence  of  computations  which  are  run  in  order  to  compute 
the  static  aerodynamic  parameters,  including  turbulent  viscous  effects, 
is  shown  in  Figure  3.  Each  block  represents  a  separate  computer  code. 
These  codes  have  been  combined  using  the  overlay  technique  on  the  BRL 
Cyber  computer.  The  two  main  codes  are  those  which  compute  three 
dimensional  turbulent  boundary  layer  development  and  three  dimensional 
inviscid  flow. 


ttauocE  or  cowuranaM 


Figure  3.  Sequence  of  Computations 


The  computation  of  the 
effects  of  viscosity  is  of  crucial 
importance  when  such  parameters  as 
roll  damping,  Magnus,  and  drag  are 
of  interest.  The  technique 
employed  here  is  a  fully  implicit, 
finite  difference  numerical  scheme 
developed  by  Dwyer ^  .  This 

technique  takes  into  consideration 
the  changes  in  direction  of  the 
cross- flow  velocity  that  occur  on 
the  side  of  the  shell  where  the 
inviscid  cross- flow  opposes  the 
surface  spin. 

The  equations  solved  are 
the  basic  equations  defining  the 
three-dimensional  compressible, 
turbulent  boundary- layer  flow  over 
a  body  of  revolution  described  by 
the  relation  r  =  r(x).  The 
coordinate  system  is  shown  in 
Figure  4. 


Figure  4.  Coordinate  System 


STUREK,  MYLIN  $  BUSH 


Continuity 


al  ^  +  37  ^  +  a?  (piJ)  =  0 


x  momentum 


-f-  3u  A3u  w  3u  w2  3rl  3  f  3u  — y 

PLU  3X  +  V37  *T  H  '  “  33d  =  -  “57  +  3y  Lu3y  -  PU  V  J 


(1) 


(2) 


4>  momentum 


-f-  3w  "3vi  w  3w  uw  3rl  1  9^e  3  f  3w  — y- J 

P[_U  37  +  v37  +  7  3*  +  -r  3ld  =  -  r  If  +  3?  \^y  '  pV  W  J  (3) 


Energy 
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PLU  a?  +  va7  +  r  afj  =  u  tx  +  F  1*  +V[\W)  \Wl  J 


pu  V 


3u 

3y 


— r-r  3w  3  U  3h  — yr -y 

-  pv  w  — —  +  —  — - pv  h 

3y  3y  Pr  3y 


(4) 


where  v  =  v  +  p'v'/p  and  the  bar  indicates  a  time- averaged  quantity. 

In  order  to  obtain  closure  of  this  system  of  equations,  the 
following  models  of  the  turbulence  terms  have  been  introduced: 

Turbulent  shear  stress 


pu  V  =  -  pv 


'■■»■§)'•(#)■«)■•  (iff 


where  e  is  introduced  as  the  turbulent  viscosity  and  the  mixing  length, 
l  =  0.09  <Stanh[ (0.4/0. 09) (y/6)] .  Van  Driest  damping  is  used  to  account 
for  the  effect  of  the  laminar  sublayer. 

Turbulent  heat  transfer 

kt  3h 

-  pv  h  =  —  y— 
cp  3y 

The  turbulent  Prandtl  number  is  introduced  as 

Pr  =  c  e/k  =  0.90 
t  p  t 
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The  three-dimensional  displacement  surface  is  not  merely  the 
vector  sum  of  the  longitudinal  and  circumferential  components  of  the 
boundary- layer  displacement. thickness.  Instead,  the  differential 
equation  derived  by  Moore 

4  iw<s;-d-5pi  *  4  toeVs3-D-8;»  ■  °  <5> 

must  be  solved  for  6*  the  three-dimensional  boundary- layer  displace¬ 
ment  thickness  where 


With  a  body  fixed  coordinate  system,  the  gas  dynamic 
equations  for  inviscid  flow  can  be  written  as 


E  +  F  +  G  +  H  =  0 
z  r  $ 


(6) 


where  the  flux  vectors  E,  F,  G,  and  H  are 


Pu, 

pv 

pw 

pv 

E  = 

pu^+  p 
puv 

F  = 

puv 

pv2+p 

G  =  - 
r 

puw 

pvw 

H  = 

puv 

p(v2-w2) 

J3UW 

_pvw 

pw2+p 

_2pvw 

These  equations  are  solved  using  MacCormack ' s ^  two-step, 
predictor- corrector  finite  difference  scheme.  The  unique  feature  of 
the  program  used  here,  which  was  developed  by  Sanders for  the  Magnus 
problem,  is  that  the  flow  field  is  computed  about  an  axisymmetric  model 
plus  displacement  surface.  Due  to  the  distortion  of  the  viscous  layer 
caused  by  interaction  of  the  surface  spin,  the  effective  aerodynamic 
shape  has  no  plane  of  symmetry. 


The  flow  field  variables  resulting  from  these  computation 
steps  have  been  developed  to  yield  the  following  aerodynamic 
coefficients--pitching  moment,  normal  force,  center  of  pressure,  Magnus 
force,  Magnus  moment,  Magnus  center  of  pressure,  form  drag,  viscous 
drag,  and  roll  damping.  The  computational  time  for  a  single  body 
configuration  and  flow  field  condition  is  approximately  ten  minutes  on 
a  CDC  7600  computer. 
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Coning  Motion  Computations 


In  order  to, compute  the  effective  pitch  damping,  the  technique 
developed  by  Schiff^8'  is  used.  This  computational  technique  relates 
the  side  force  on  a  body  undergoing  a  coning  motion  about  the  CG 
location  to  the  pitch  damping  (CM  +  CM> ),  see  Figure  5. 

q  a 


The  numerical  technique 
is  MacCormack's16^  predictor- 
corrector,  explicit  marching  scheme. 
This  computation  involves  the 
solution  of  the  Euler  equations 
including  terms  for  Coriolis 
(2pxfjxv)  and  centrifugal  [pSix(jixf)] 
forces  in  a  body  fixed  coordinate 
system.  For  this  case,  the  H  vector 
in  equation  6  becomes 


Figure  5.  Coning  Motion  About  Center  of  Gravity 


PV 


puv+pr[2(u)2W-w3v)+wiui2r-x(u)£+u)§)  ] 

p  (v2-w2)+pr[2(w3u-u>iw)+u>iu>2X-r(u>f +u>3)  ] 


_2pvw+pr[2(u>iv-u)2u)+u>3  (u>2r+u>is)] 


where  u>i,  (02,  and  <*>3  are  the  components  of  the  angular  velocity  vector 
6  resolved  in  the  z,  r,  and  <j>  directions,  respectively. 


For  the  case  of  a  steady  coning  motion,  the  flow  field  is 
time-invariant  in  the  body- fixed  coordinate  system.  The  effective 
pitch  damping  (C^  +  C^#)  is  determined  using  the  relation 

q  a 


* 
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C  *  sina(CM  ♦ 
9  q 


(7) 


where  C  =  side  force  at  coning  rate  0  and  effective  angle  of  attack 

Ye 

o,  which  is  valid  for  small  values  of  a  and  6.  Thus  a  dynamic 
aerodynamic  parameter  is  determined  using  a  steady  flow  field 
computation.  This  is  a  potentially  very  useful  tool  for  the  exterior 
ballistician.  The  computation  time  is  approximately  90  seconds  on  a 
CDC  7600  computer  for  the  body  configurations  in  this  study. 


III.  RESULTS 


Comparisons  to  Experiment 

Detailed  comparisons  of  the  computations  to  experimental  data 
for  turbulent  boundary  layer  profile  characteristics,  wall  pressure 
measurements  and  Magnus  force  are  reported  in  Ref.  (9).  Comparisons 
shown  here  will  be  limited  to  the  aerodynamic  coefficients  of  interest. 

Examples  of  comparisons  of  the  computed  results  to  experimental 
data  are  presented  in  Figures  6  through  10.  The  comparisons  for  pitch 
plane  static  parameters  shown  in  Figures  6,  7  and  8  indicate  excellent 
agreement.  The  limited  comparison  for  Magnus  in  Figure  9  indicates 
acceptable  agreement  if  allowance  is  made  for  the  small  magnitude  of 
the  Magnus  effect  and  the  variance  between  the  wind  tunnel  and  range 
experimental  measurements.  A  comparison  between  computation  and 
experiment  for  pitch  damping  is  shown  in  Figure  10.  The  experimental 
point,  which  is  for  an  L/D  of  5.12  and  cone  angle  of  9.52°,  shows 
excellent  agreement  with  the  trend  of  the  computed  results.  In  general, 
it  is  felt  that  the  numerical  computations  do  provide  an  accuracy  for 
the  aerodynamic  coefficients  that  is  within  the  uncertainty  of  our 
ability  to  determine  these  coefficients  experimentally.  However,  it  is 
felt  that  a  broader  scope  of  comparison  for  the  aerodynamic  coefficients 
between  experiment  and  computation  is  of  interest  and  increased  effort 
to  accomplish  this  is  underway. 
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Figure  6.  Pitching  Moment,  Comparison  with  N*™  7-  ^enter  of  Pressure.  Comparison  with 
Experiment  Experiment 


Figure  8.  Normal  Force,  Comparison  with 
Experiment 


Figure  9.  Magnus  Moment,  Comparison  with 
Experiment 
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Figure  10.  Pitch  Damping,  Comparison  with 
Experiment 
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Parametric  Comparisons 

Examples  illustrating  the  parametric  results  are  shown  in 
Figures  11  through  26.  The  series  of  comparisons  shown  in  Figures  11 
through  20  illustrates  an  example  for  each  aerodynamic  coefficient 
computed  in  this  study.  The  case  chosen  is  the  SOC  model  for  a  total 
length  of  six  calibers  and  for  ogive  lengths  of  two,  three,  and  four 
calibers.  The  aerodynamic  coefficients  are  plotted  versus  Mach  number 
for  atmospheric  free  stream  launch  conditions  assuming  an  adiabatic 
wall  temperature  boundary  condition.  These  comparisons  show,  for  a 
fixed  body  length,  that  configurations  with  long  slender  ogives  have 
reduced  pitch  damping,  less  drag,  and  a  reduced  Magnus  moment  compared 
to  bodies  with  shorter  ogive  lengths. 


Figure  13.  Normal  Force,  Parametric 
Comparison,  SOC,  L/D  *  6 
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Figure  14.  Magnus  Moment,  Parametric 
Comparison,  SOC,  L/D  =■  6 
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Examples  are  shown  in  Figures  21  through  23  illustrating  the  effects  of 
variations  in  ogive  shape  for  fixed  forebody  and  total  projectile 
lengths.  These  comparisons  show  that  pitching  moment,  Magnus  moment, 
and  pitch  damping  are  increased  as  ogive  bluntness  is  increased. 


I/O  •  6 ,  ln  *  3 
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The  final  sequence  of  parametric  comparisons  is  shown  in  Figures  24 
through  26  where  the  effect  of  varying  the  body  length  is  shown  for 
fixed  ogive  shape.  These  figures  show  that  pitching  moment,  Magnus 
moment,  and  pitch  damping  are  all  increased  as  the  body  length  is 
increased. 


SOC,  L N  *  3 

CS. Ar  061  BEHIND  NOSE 


SOC  ,  IN  .3 
CS.AT  0.61  BEHIND  NOSE 


Figure  24.  Pitching  Moment,  Parametric 

Comparison,  SOC,  L..  =  3 
N 


Figure  25.  Magnus  Moment,  Parametri 
Comparison,  SOC,  L..  =  3 
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Figure  2b.  Pitch  Damping,  Parametric 
Comparison,  SOC,  L^.  =  3 


STORE K,  MYLIN  $  BUSH 


The  comparisons  shown  represent  a  small  fraction  of  the 
potential  comparisons  possible  from  the  total  data  base  generated.  The 
intent  here  has  been  to  illustrate  the  capability  of  the  computation 
techniques  rather  than  develop  any  conclusion  as  to  the  relative 
superiority  of  any  particular  configuration.  This  study  is  part  of  a 
continuing  effort  that  is  being  expanded  to  include  boattail 
configurations  and  a  wider  Mach  number  range- -transonic  velocities  are 
of  particular  interest. 

IV .  SUMMARY 

A  computational  aerodynamics  parametric  study  has  been 
described  in  which  advanced  numerical  techniques  for  computing  three- 
dimensional  inviscid  and  turbulent  viscous  supersonic  flow  fields  have 
been  used.  A  comprehensive  data  base  has  been  generated  for  cone- 
cylinder,  tangent-ogive-cylinder,  and  secant-ogive-cylinder 
configurations.  Comparisons  between  the  computed  results  and  experiment 
have  provided  verification  of  the  computational  techniques.  Comparisons 
of  the  computed  results  for  differing  body  configurations  have 
established  the  ability  of  the  computation  techniques  to  distinquish 
the  effects  of  body  configuration  on  the  aerodynamic  coefficients. 

The  comparisons  for  this  study  represent  the  first 
comprehensive  use  of  advanced  flow  field  computation  techniques  to 
develop  a  parametric  aerodynamic  analysis  which  includes  turbulent 
viscous  effects.  Of  particular  interest  are  the  computations  of  Magnus 
effects,  which  are  accomplished  in  a  conceptually  exact  manner,  and  the 
computations  of  pitch  damping.  Although  the  technique  of  computing 
pitch  damping  from  a  steady  coning  motion  has  been  available  for  several 
years,  this  study  is  the  first  known  comprehensive  application  of  this 
technique  to  ogive-cylinder  bodies. 
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NOMENCLATURE 

A  =  reference  area  =  ttD2/4 

CDBL  =  viscous  drag  =  (// t  cose_dS)/qA 

CDF  =  total  drag  =  (//pwsin0gdS)/qA  +  CDBL 

CLP  =  roll  damping  =  (J/rx^dS)/ (qAD  PD/V) 

CM  =  pitching  moment  =  (//z p  coS(|>cos0DdS)/qAD 

CMA  =  pitching  moment  coefficient  =  CM/oc 

CMPA  =  Magnus  moment  coefficient  =  CMY/(PD/V/a) 

CMY  =  Magnus  moment  =  [//( zp  sin<f>cos0D  +  zt  ,  cos^cosO,,  +  zApsin<|>cos6 

W  D  <p  D 

+  zTxsin<|>sin0g)dS]/ (qAD) 

CN,  C^  =  normal  force  =  (J7pwcos<t>cos6gdS)/qA 

O 

CMQ  +  CMA=  pitch  damping  =  Cn*/sino 

CNY,CV  =  Magnus  force  =  f//(p  sin<j>cos0D  +  t , cos<j>cos0_  +  Apsin<j>cos0„  + 

I  W  D  <p  D  D 

t^s  in<(>  s  in0  B)  dS  ] /qA 
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i 


CYe 

CPN 

CPY 

D 

P 


u,  v,  w 
V 


x 

y.  y 

z 

Ap 

a 


6 


•  i 

=  side  force  in  coning  motion  =  (//pwsin<f>cos0gdS)/ (qA0)  i 

=  center  of  pressure  =  CM/CN 

=  Magnus  center  of  pressure  =  CMY/CNY  ] 

=  diameter  of  model 

=  spin  rate,  rad/s  j 

=  free  stream  dynamic  pressure  =  (pooV^)/2  j 

=  local  radius  of  model  i 

=  Reynolds  number  based  on  model  length 

=  surface  area 

=  velocities  in  boundary- layer  coordinates 
=  velocity  along  model  trajectory 
=  surface  coordinate  in  longitudinal  direction 
=  coordinate  perpendicular  to  local  surface 
=  cylindrical  coordinate  along  model  axis 
=  centrifugal  pressure  gradient  contribution  to  side  force 
=  effective  angle  of  attack  for  coning  motion 
=  longitudinal  velocity  wall  shear 

=  circumferential  velocity  wall  shear 

=  local  slope  of  body  surface  5 

=  coning  rate 


